MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURtAU  OF  STANDARDS- 19S3-A 


.  . . * -..■■{W.MAORBIIIWiWIW— 


l 


T 


UNCLASSIFIED _ . 

igfaftTTVTLASSlftdATHDN  6F  Vhi4  PA5F 


ADAIS2.SSF 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 
4.  PERFORMING  ORGANIZATION  REPORT  NUMBERS) 


6-1.  NAME  OF  PERFORMING  ORGANIZATION 

UNIVERSITY  OF  ILLINOIS 


/LOOMISCIJ^BO^TORY  OF^HYSICS 
I110W  GREEN  STREET 
URBANA  IL  61801 


6b.  OFFICE  SYMBOL 
(If  applicable) 


8a.  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 
same  as  7a 


8b.  OFFICE  SYMBOL 

(If  applicable) 

_n£ _ 


lb.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 
Approved  for  Public  Release; 
Distribution  Unlimited 


S. 


7a.  NAME  OF  MONITORING  ORGANIZATION 

AFOSR/NP 


7b  tnd  ZIPCo*> 

BOLLING  AFB  DC  20332-6448 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

AFOSR  84-0384 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

ACCESSION  NO 

same  as  7b 

61102F 

2301 

A8 

1 1 .  TITLE  (Include  Security  Cassifkation) 

Thermodynamics  of  Strain-confined  Paraexcitons  in  CU(2)0  (U) 


12.  PERSONAL  AUTHOR(S) 

D  P  Trauernicht.  J  P  Wolfe,  &  A  Mysyrowicz 


OF  REPORT 


t3,b  J^p^Q  28  Feb  ^  OATE  gFbREgyT  (Year,  Month,  Day)  jlS  ^.E  COUNT 


16.  SUPPLEMENTARY  NOTATION 


1 17.  COSATI  CODES  1 

FIELD 

GROUP 

SUB-GROUP 

18.  SUBJECT  TERMS  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 
Paraexcitons,  Phonons,  B-E  Condensation 


fteA“SESy8^  confined  to  a  parabolic  potential  well 

was  examined.  From  time  and  space  resolved  studies  of  the  ortho  and  para  luminescence, 
plausible  causes  for  the  saturation  of  paraexciton  density  in  the  strain  well  were  deduced. 


20.  DISTRIBUTION  /  AVAILABILITY  OF  ABSTRACT 
OUNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT.  □  OTIC  USERS 


21 .  ABSTRACT  SECURITY  CLASSIFICATION 
UNCLASSIFIED 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

PI  Robert  J. *  Bdrlter _ 


[22b.  TELEPHONE  (Include  Area  Code) 

202/767-5011 


22c.  OFFICE  SYMBOL 
NP 


DO  FORM  1473, 84  MAR 


83  APR  edition  may  be  used  until  exhausted. 
All  other  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 

87  7 


1  126 


Thermodynamics  of  Strain-confined  Paraexcitons  in  Ci^O- 


D.  P.  Trauernlcht  and  J.  P.  Wolfe 
Department  of  Physics 

University  of  Illinois  at  Urbana-Champaign 
1110  W.  Green  Street 
Urbana,  IL  61801 


and 


A.  Mysyrowicz 

Group  de  Physique  des  Solides 
Ecole  Norma le  Superieure 
Universite  Paris  VII 
Tour  23,  F-75005 
Paris,  France 


Abstract 


Accesion  For  j 

NTIS  CRA&I 

DTIC  TAB 

□ 

Unannounced 

□ 

Justification 

By 

Diiti  ib'jtion/ 

Availability  Codes 


Dist 


IM 


Avail  j  djor 
Special 


The  thermodynamic  behavior  of  long^iived  paraexcitons  confined  to  a 
parabolic  potential  well  is  examined.  The  potential  well  is  produced  by  a 
Hertzian  contact  stress.  A  wavelength^cunable  dye  laser  is  used  to  create 
excitons  directly  in  the  potential  well  or  at  any  other  localized  point  inside 
the  crystal.  Spectral  and  spatial  distributions  of  the  exclton  recombination 
luminescence  are  measured  for  CW  and  pulsed  excitation.  The  possibility  of 
Bose^lnstein  condensation  (BEC)  of  these  long-lived  excitons  is  examined  both 
theoretically  and  experimentally.  We  calculate  the  spectral  and  spatial 


distribution  of  luminescence  from  a  gas  of  non-interacting  particles  in  a 


three£dlnenslonal  harnonlc^&sclllator  veil#  The  results  are  narkedly 

c  °  - 

different  for  direct  (noQphonon)  and  indirect  (phononf&ssisted) 

recombination.  The  calculated  spectra  are  compared  to  the  data  for  moderate 
CW  excitation^^t^T^^  *2-4.2  K^js,Taken  alone,  the  mx^l&onon  spectra 
suggest  that  th*  excltonic  gas  is  In  the  quantum  regime;  however,  this 
conclusion  is  shown  to  be  inconsistent  with  the  estimated  density  of  the 


2 


l 


A 

gas.  A  consistent  Interpretation  of  all  spectral  and  spatial  distributions  is 
possible,  however,  if  one  assumes  Maxwell-Boltzaann  statistics  and  takes  into 
account  the  rapidly  changing  paraexclton  intensity  with  applied  stress.  From 
time  and  space -resolved  studies  of  the  ortho  and  para  luminescence,  plausible 
causes  for  the  saturation  of  paraexclton  density  in  the  strain  well  are 
deduced.  First,  an  Auger-like  recombination  of  colliding  paraexcltons  seems 
to  limit  their  density  —  an  idea  supported  by  the  observed  power  dependence 
of  ortho  and  paraexclton  signals.  Secondly,  an  anomalously  slow 
thermalizatlon  of  strain-confined  paraexcltons  is  observed. 


PACS  Nos:  71.35.-fz 
78.55.-m 
Typed  by:  T.  Watts 
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I.  Introduction 

The  possibility  of  Bose-Elnstein  Condensation  (BEC)  of  excitons  in 
semiconductors  has  been  actively  pursued  for  many  years.  The  exdton,  or 
bound  electron-hole  pair,  is  the  basic  electronic  excited  state  of  an 
intrinsic  semiconductor.  Its  light  effective  mass  is  a  favorable  factor  for 
the  occurrence  of  BEC  at  moderate  particle  densities  and  relatively  high 
temperatures.  Specifically,  for  an  ideal  gas  of  particles,  a  macroscopic 
occupation  of  the  ground  state  occurs  at  densities  above  nc  -  (6.2  x  10 ^cm-^) 
g(m*/m  )T^^,  which  yields  nc  ■  2  x  lQ^crn*^  at  T  ■  10  K  for  an  exciton  mass 
m*  equal  to  the  free  electron  mass  mQ  and  a  degeneracy  factor  of  g  >  1.  The 
observation  of  this  quantum  statistical  effect,  however,  has  proven  quite 
elusive,  mainly  due  to  coupling  effects  from  exciton-exciton  interactions. 

For  example,  excitons  may  undergo  a  Mott  transition  to  an  electron-hole  plasma 
or  condense  into  electron-hole  droplets1  before  nc  is  reached.  Recent  studies 
have  therefore  concentrated  on  crystals  for  which  the  Mott  and/or  spatial 
condensation  densities  are  expected  to  be  higher  than  nc. 

One  excitonic  system  which  was  reported  to  have  shown  quantum  statistical 

o 

effects  is  a  gas  of  blexcltons  in  CuCl.  Chase  et  al.  used  two-photon 

absorption  of  a  high  Intensity  laser  to  create  a  large  initial  density  of 

blexcltons  at  twice  the  photon  wavevector.  When  a  weak  density  of  probe 

particles  was  added  to  the  system,  a  redistribution  of  these  probe  particles 

in  k-space  strongly  suggested  the  presence  of  a  Bose  condensate.  Also, 

excitons  in  Ge  have  been  reported  to  show  quantum  statistical  effects 

characteristic  of  Bose  particles.  For  unlaxially  compressed  Ge  in  a  magnetic 

•» 

field,  Timofeev  et  al.  observed  a  narrowing  of  the  exciton  luminescence 
llneshape  with  increased  excitation,  which  they  interpreted  as  evidence  for 
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Bose-Einstein  statistics.  It  was  suggested  that  condensation  might  occur  in 
this  system  at  temperatures  below  T  -  1  K. 

Another  crystal  which  has  been  suggested  as  a  good  candidate  for  the 
observation  of  excitonic  BEC^  is  C^O.  In  this  material  the  lowest  exciton  level  (n 
-  1)  is  split  into  two  states:  a  triply-degenerate  orthoexciton  state  and  a 
nondegenerate  paraexciton  state  lying  at  12  meV  lower  energy.  Luminescence 
experiments  have  been  reported  by  Hulin  et  al. ^  showing  evidence  of  Bose  statistics 
for  the  orthoexcitons.  As  the  excitation  level  was  increased,  both  the  quantum 
degeneracy  and  the  temperature  of  the  excitonic  gas  increased.  Hulin  et  al. 
concluded  that  the  critical  density  nc  for  BEC  was  reached  at  their  highest 
excitation  levels  and  T  ■  20  K.  A  subsequent  theoretical  study  was  carried  out  by 
Haug  and  Kranz,**  proposing  that  the  data  of  Ref.  5  contained  evidence  that  the 
particle  density  was  in  fact  exceeding  nc«  They  noted  the  presence  of  a  low  energy 
tall  in  the  data  which  could  be  attributed  to  the  scattering  of  excitons  from  a 
condensed  fraction  of  particles. 

In  this  paper,  we  present  a  study  of  excitons  in  Cu20  confined  in  a  stress- 

2 

induced  parabolic  potential  well,  with  potential  energy  V(r)  -  or  .  We  employ  a 
Hertzian  stress  configuration  to  confine  the  excitons  to  a  controlled  volume  Inside 
the  crystal.  This  geometry  looks  appealing  for  the  achievement  of  BEC  of  para- 
exdtons,  since  in  principle  one  should  be  able  to  collect  a  high  density  of  cold 
particles,  away  from  the  hotter  excitation  spot  and  from  undesirable  surface 
effects.  We  find  that  the  paraexcltons  exhibit  a  quasi-equilibrum  in  the  parabolic 
potential  well.  This  is  possible  due  to  the  long  lifetlma  of  paraexcltons,  owing  to 
the  fact  that  the  radiative  recombination  of  excitons  is  forbidden  in  Cu20,  despite 
the  direct  band  gap  structure  of  this  material.  It  was  previously  reported  that 
naturally  grown  crystals  of  Cu20  exhibited  paraexciton  lifetimes  exceeding  a 
microsecond.7  We  employ  samples  of  the  same  material  for  our  studies. 


.J 
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Generally  Bose-Elnsteln  condensation  la  described  for  an  Ideal  gas  of 
free  particles  In  a  region  of  constant  potential.  Significant  changes  are 
expected,  however,  for  Bose  particles  In  an  external  potential.  Therefore  we 
begin  by  examining  theoretically  the  quantum-statistical  behavior  of  bosons  In 
a  parabolic  potential  well.  Predictions  are  made  concerning  both  the  spatial 
and  energy  distribution  of  the  particles.  Unlike  the  usual  BEC  for  a  constant 
potential— which  is  is  a  k- space  condensation — BEC  In  a  harmonic-oscillator 
potential  Implies  a  real-space  condensation.  We  show  that  the  spatial  and 
spectral  luminescence  profiles  are  expected  to  be  quite  different  for  no- 

fi 

phonon  and  phonon-assisted  excitonic  recombination  lines.  The  spatial 
distribution  of  phonon-assisted  luminescence  reflects  the  actual  particle 
density. 

As  will  be  shown,  interpretation  of  the  strain-confined  paraexciton 
luminescence  spectra  in  terms  of  the  above  theory  presents  some  ambiguities. 
The  measured  spectra  are  found  to  be  in  good  agreement  with  the  theoretical 
predictions  and  could  be  interpreted  as  due  to  an  ideal  gas  of  particles  in 
the  quantum  regime,  close  to  BEC.  However,  calibration  of  the  experimental 
particle  densities  in  the  well  precludes  this  possibility.  This  discrepancy 
may  be  resolved  by  taking  into  account  the  strong  stress  dependence  of 
paraexciton  luminescence  intensity  and  using  classical  statistics. 

In  attempting  to  achieve  high  densities  of  paraexcitons ,  we  encountered 
some  limitations  associated  with  exciton  -  exclton  collisions.  Under 
increasing  laser  excitation,  the  paraexciton  signal  was  found  to  behave  in  a 
very  aubllnear  fashion.  This  subllnear  behavior  is  explained  in  terms  of  an 
Auger-type  recombination  channel,  where  two  excltons  collide  resulting  in  the 
lose  of  one  electron-hole  pair  and  the  creation  of  a  free  electron  and  a  free 
hole  which  subsequently  thernellae  to  form  an  orthoexciton.  The  final  section 
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of  this  paper  characterizes  these  thermallzation  processes  by  tine-resolved 
luminescence  measurements . 

II.  Bose  Statistics  in  a  Parabolic  Well 

We  wish  to  predict  the  luminescence  spectrum  and  spatial  distribution  of 

excitons  that  are  confined  in  a  3-dimensional  harmonic  oscillator  potential, 

V(r)  *  ar  .  If  it  is  assumed  that  they  behave  as  a  classical  gas  and  that  the 

matrix  element  for  radiative  recombination  is  Independent  of  wavevector  and 

spatial  position,  the  spatial  profile  can  be  obtained  quite  simply.  For  a 

classical  gas,  the  internal  chemical  potential  is  pint  =  kgT  ln(n/nq),  where  n 

■  n(r)  is  the  density  and  n  *  (m  k_T/2wH  )  is  the  quantum  density.  The 

4  a 

total  chemical  potential,  p,  is  the  sum  of  the  internal  and  external  chemical 
potential  (Pext  ■  or2).  Since  u  «*  constant  in  the  potential  well,  one  obtains 

kgT  lntnCoj/n^)  -  kgT  ln(n(r)/nq)  +  or2  (1) 

which  has  the  solution 

n(r)  *  n(o)  exp(-ar^/k  T).  (2) 

D 

One  way  of  obtaining  the  luminescence  lineshape  under  the  above 
assumptions  is  to  note  that  a  harmonic  oscillator  has  energy  levels  Es  **  (s  + 
3/2)Ku  with  u  ■  (2a/m  )  '  and  s  -  Each  level  has  a  degeneracy  of  (s 


+  l)(s  +  2)/2.  Since  the  energy  levels  are  equally  spaced,  the  density  of 
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states  is  given  by  the  degeneracy.  For  typical  experiments  tiu  »  10  eV  « 

kgT,  implying  that  most  particles  are  in  states  s  »  1.  The  density  of  states 
2 

then  varies  as  E  .  Assuming  a  Maxwell -Boltzmann  population  of  these  states, 
the  luminescence  lineshape  is  predicted  to  be, 

1(E)  «  E2  exp(-E/kBT  ).  (3) 

As  the  exciton  gas  approaches  degeneracy  (n  ■+■  n^),  the  above  equations  do 
not  hold.  Indeed,  there  is  not  a  closed-form  expression  for  the  chemical 
potential  of  degenerate  Bose  gas  as  a  function  of  density.  One  method  of 
calculating  the  spatial  profile  is  to  sum  over  the  probability  distribution 
for  harmonic  oscillator  wavefunctions ,  |<j>|  ,  weighted  by  the  probability  of 
occupation  of  the  state,  (exp[(E  -  u)/k_T]  -  1}  1  .  Alternatively,  one  can 

8  a 

consider  the  confined  gas  as  a  sum  of  small  volumes  with  nearly  uniform 
density,  which  is  the  approach  we  now  take.  The  results  depend  upon  whether 
the  observed  recombination  light  is  due  to  a  no-phonon  or  phonon-assisted 
transition.  The  no-phonon  transition  is  constrained  (locally)  to  the  K  »  0 
selection  rule,  whereas  the  phonon-assisted  transition  can  sample  a  full  range 
of  exciton  wavevectors. 


Phonon-assisted  transitions 

Suppose  that  the  spatial  profiles  and  luminescence  lineshapes  are  given 
by  the  sum  of  local  lineshapes,  which  are  dependent  on  the  spatial  position. 
For  phonon-assisted  recombination  in  which  the  matrix  element  is  assumed 
Independent  of  wavevector,  the  local  lineshape  is  given  by 


(4) 


rr  ^  Aoc<^> 

I  (E,r)  «  - 

exp[(Eioc  -  Mloc)/V]  -  1 

where  ^loc(E,r)  ■  (E  -  is  the  local  density  of  states,  Ejqc  “ 

2  2 
E  -  or  is  the  local  kinetic  energy,  and  P^oc  “  W  “  ar  *s  local  internal 

chemical  potential.  The  total  chemical  potential  p  is  independent  of  position 

r.  The  spatial  profile  is  obtained  by  integrating  Eq.  4  over  energy.  Note 

that  E  must  be  greater  than  ar  •  This  yields 


I(r)  « 


00 

l 

n=l 


exp[-n(ar2  -  pJ/kgT] 
__ 

n 


(5) 


If  the  system  is  not  degenerate  (p  is  large  and  negative),  this  reduces  to  the 
form  of  Eq.  2,  because  only  the  first  term  in  the  sum  will  contribute. 
Similarly,  the  luminescence  lineshape  is  obtained  by  integrating  Eq.  4  over 
all  space.  This  yields 


1(E)  «  - - -  ,  (6) 

exp  f (E  -  p)/kgT]  -  I 


which  reduces  to  the  fora  of  Eq.  3  for  a  nondegenerate  system. 

Shown  in  Fig.  la  is  Eq.  6  for  various  values  of  p  at  T  *  2K.  In  Fig.  lb 
is  Eq.  5  integrated  over  the  yz  plane  ("slit  scans")  for  the  same  values  of  p 
and  T  as  above.  These  calculations  display  very  little  change  in  the  spectral 
lineshape  or  the  spatial  profile  when  degeneracy  is  assumed.  The  width  of  the 
spectral  lineshape  does  not  change  noticeably,  but  the  peak  of  the 


distribution  does  shift  somewhat  to  lower  energy.  The  spatial  profile  narrows 
only  slightly. 

Because  the  ground  state  is  excluded  in  the  above  calculations,  a 
schematic  representation  of  its  effect  is  shown  in  the  bottom  traces  of  Figs, 
la  and  lb.  The  ground  state  population  is  simply  included  as  a  spike  on  the 
low  energy  edge  of  the  spectrum  and  in  the  center  of  the  parabolic  well, 
respectively.  This  illustrates  an  important  aspect  of  strain-confined 
excltons:  BEC  into  the  lowest  energy  level  (the  harmonic  oscillator  ground 
state)  implies  a  spatial  condensation.  This  is  unlike  the  constant-potential 
case,  for  which  there  is  a  condensation  to  K  =  0  in  wavevector  space,  not  a 
spatial  condensation. 

A  proper  calculation  of  the  particle  distribution  with  the  inclusion  of 
the  ground  state  needs  to  be  done.  Specifically,  the  low  energy  spike  on  the 
spectral  lineshape  could  be  broadened  due  to  interactions  in  the  condensate. 
The  recombination  event  may  involve  simultaneous  emission  or  absorption  of  a 
coherent  ground  state  excitation,  as  pointed  out  by  Haug  and  Kranz.^ 

The  above  analysis  indicates  that  for  phonon-assisted  transitions,  which 
sample  the  entire  exciton  energy  distribution,  the  approach  to  quantum 
statistics  will  not  be  easily  seen.  Only  when  the  system  reaches  the  BEC 
density  and  there  is  a  partial  condensate  will  there  be  definite  signs  of 
degeneracy. 


Direct  transitions 

In  the  case  of  a  direct-gap  semiconductor,  no-phonon  transitions  are  also 
observable.  No-phonon  transitions  sample  the  states  near  K  -  0,  specifically 
those  with  wavevectors  equal  to  the  emitted  photon  wavevector.  The  photon 
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momentum  Is  usually  small  compared  to  the  spread  of  momentum  in  the  system  of 
particles.  One  can  calculate  the  luminescence  distributions  by  summing  local 
lineshapes  again.  Assuming  that  the  no-phonon  transition  is  vertical  (i.e. 
ignoring  the  photon  momentum),  the  local  lineshape  is  given  by 


Il.e(E'r)  - 


S<Eloc’ 


expKEloc  ■  Mloc)/kBT1  '  1 


(7) 


where  the  6-function  picks  out  the  local  zero  kinetic  energy  state,  with  the 
number  of  particles  in  this  state  being  given  by  the  Bose  factor.  To  obtain 
the  spatial  profile,  one  integrates  Eq.  7  with  respect  to  energy.  This  is 
quite  easily  done  because  of  the  6 -function,  and  one  obtains 


I(r)  « 


_ 1 _ 

exp[(ar2  -  yJ/kgT]  -  I 


(8) 


Eq.  8  reduces  to  the  form  of  Eq.  2  for  a  nondegenerate  system  (y  large  and 
negative).  The  luminescence  lineshape  is  obtained  by  integrating  Eq.  7  over 
all  space,  which  yields 


1(E)  « 


E1/2 

'  • 

exp{(E  -  y)/kgT]  -  1 


(9) 


1/2 

The  lineshape  is  significantly  different  from  Eq.  6  in  that  Eq.  9  has  an  E 

2 

multiplier  while  Eq.  6  has  an  E  multiplier. 
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Shown  In  Fig.  2a  is  Eq.  9  for  various  values  of  y  at  T  ■  2K.  In  Fig.  2b 
is  Eq.  8  integrated  over  the  yz  plane  for  the  same  values  of  p  and  T.  In  this 
case  of  no-phonon  transitions,  there  is  a  significant  change  in  the 
luminescence  lineshape  and  spatial  profile  as  the  system  becomes  degenerate. 
This  is  due  to  the  local  K  =  0  selection  along  with  the  fact  the  density  is 
higher  in  the  center  of  the  well,  implying  higher  degeneracy  there.  A  more 
degenerate  local  lineshape  (remember  Uioc  =  U  ~  ar  )  has  a  larger  number  of 
particles  in  the  E^oc  “  0  state.  This  effect  implies  a  sharper  measured 
spatial  profile  than  for  the  phonon-assisted  spatial  distribution,  which 
samples  all  K-vectors.  Shown  as  the  bottom  trace  in  each  part  of  Fig.  2  is 
schematically  what  would  be  expected  for  a  partially  condensed  system.  Both 
Eqs.  8  and  9  contain  singularities  when  p  *  0.  These  are  indicated  by  the 
arrows  pointing  upward. 

From  this  analysis  of  the  no-phonon  lineshape  and  spatial  profiles,  it 
appears  that  the  no-phonon  line  should  be  useful  when  searching  for  Bose 
statistical  effects  for  excitons  confined  in  a  parabolic  well.  There  are 
gradual  changes  in  the  lineshape  and  spatial  profile  which  signal  the  onset  of 
quantum  degeneracy,  in  contrast  to  the  phonon-assisted  transitions.  However, 
the  phonon-assisted  transition  may  provide  a  more  definitive  identification  of 
condensation. 


Contrasts  with  constant-potential  case 

For  BEC  in  a  harmonic  oscillator  well,  the  parameter  which  determines  the 
occurence  of  condensation  at  a  given  temperature  is  the  total  number  of 
particles  in  the  well  rather  than  the  density.  This  is  seen  by  calculating 
the  total  number  of  particles  in  the  well 
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(s  +  l)(s  +  2 )/2 

N  -  g  l  -  (10) 

s  exp{  [  (n  +  3/2)Kui  -  p]/kgT}  -  1 

where  g  Is  an  additional  degeneracy  factor  (e.g.,  a  spin-like  degeneracy). 
Setting  |i  »  0  to  obtain  the  maximum  number  that  can  be  held  in  the  excited 
states  in  equilibrium,  assuming  s  »  1  as  before,  and  changing  the  sum  over  a 
to  an  integral  over  ds,  one  obtains 


N  -  g  (1.202)  (k _T/fo>)3 

C  B 


(ID 


This  critical  number  can  be  "converted"  to  a  critical  density  by  dividing  by 

the  volume  of  the  well  occupied  by  excitons,  as  determined  by  the  full-width 

at  half-maximum  (FWHM)  of  the  spatial  distribution;  i.e., 

nWeH  *  N  /fl,  with  the  effective  volume  given  by 
c  c 


fl  =  (J)  (2.77kgT/a)3/2  . 


(12) 


The  factor  2.77  comes  from  choosing  the  FWHM  of  the  density  profile  rather 
than,  say,  the  1/e  point.  The  simple  Gaussian  spatial  profile  result  is 
assumed  because  there  is  only  a  minor  change  in  the  width  of  the  spatial 
distribution  of  the  particles  in  the  excited  states  when  approaching  Bose 
degeneracy,  as  shown  above.  (The  phonon-assisted  transition,  Fig.  lb,  mimics 
the  local  density  whereas  the  no-phonon  transition.  Fig.  2b,  does  not.)  This 
results  in  a  critical  density  given  by 


well 

n 


(T) 


0.176 


<V> 


3/2 


(13) 
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Note  that  the  well  parameter  a  does  not  appear  In  Eq.13.  This  may  be  compared 

Q 

with  the  result  for  "particles  In  a  box",  which  Is 

3 /2 

nfree  (T)  -  0.166  —S~ -  (k  T)3/2  .  (14) 

c  K3  8 

There  Is  little  difference  between  the  densities  In  Eqs.  13  and  14.  For  this 
reason,  one  can  still  speak  of  a  critical  density  for  BEC  of  excitons  In  a 
parabolic  well,  with  this  density  being  the  same  as  that  for  "free"  particles. 


III.  Strain  Confinement  of  Paraexcltons 

To  produce  a  parabolic  well  inside  the  crystal,  an  inhomogeneous  stress 
is  applied  to  the  sample  by  pressing  a  rounded  plunger  against  one  face  of  the 
sample.  This  Hertzian  stress  geometry  has  been  described  elsewhere  for  the 
indirect-gap  semiconductors  Si  and  Ge.*  ’ ^  A  maximum  of  shear  stress  is 
produced  at  a  finite  distance  into  the  sample.  As  discussed  below,  the 
excitons  are  attracted  to  this  shear-stress  maximum  (SSM)  and  confined 
there.  The  parabolic  dependence  of  the  external  potential  acting  on  the 
excitons  comes  about  from  a  Taylor  expansion  of  the  stress  distribution  near 
the  SSM. 

Excitons  are  attracted  to  a  SSM  if  their  energy  is  lowered  by  the 
application  of  a  shearing  stress.  For  paraexcltons  in  (h^O,  the  application 
of  an  uniaxial  stress  along  a  four-fold  axis  lowers  their  energy  as  reported 

12  13 

in  an  earlier  paper.  Waters  et  al.  have  shown  through  second  order 
perturbation  theory  that  the  paraexciton  energy  as  a  function  of  stress  can  be 
expressed  in  the  following  form: 
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Epara  «  A  +  Ba  -  Ca2  (15) 

where  a  Is  Che  magnitude  of  the  uniaxial  stress.  Each  term  is  easily 

identified.  A  is  the  unperturbed  paraexciton  energy.  Bo  is  due  to  the 

o 

dilational  term  in  the  strain  Hamiltonian,  and  Co  is  from  a  term  of  lowered 

2 

symmetry  in  the  strain  Hamiltonian.  It  is  the  dominance  of  the  Co  term  at 
high  applied  stresses  that  makes  the  SSM  an  energy  minimum  for  the 
paraexcitons  because  this  term  incorporates  non-dilational  (shear)  strains. 

Shown  in  Fig.  3  are  two  spectra.  The  lower  trace  is  a  photoluminescence 
spectrum  from  a  crystal  with  no  stress  applied  and  excited  by  an  argon-ion 
laser.  Xq  denotes  orthoexciton  no-phonon  line;  the  other  lines  are  optical 
phonon  replicas  of  the  orthoexciton  line  and  are  identified  by  the  symmetry  of 
the  phonon  emitted.  One  exception  is  a  phonon  replica  of  the  paraexciton 
labeled  Xp  -  The  upper  trace  is  a  spectrum  taken  from  the  SSM  at  an 

approximate  stress  of  3.6  kbar.  Note  that  the  wavelength  scale  is  shifted. 

An  additional  line  appears  under  the  application  of  a  symmetry-lowering  stress 
and  is  identified  as  the  no-phonon  line  of  the  paraexciton,^  Xp.  This 
spectrum  is  obtained  by  resonantly  pumping  (with  a  tunable  dye  laser)  the 
down-shifted  phonon-assisted  orthoexciton  absorption  band  at  a  spatial 
position  "below"  the  well.  The  longer-lived  paraexcitons  drift  into  the 
well.  The  appearance  of  a  relatively  small  orthoexciton  signal,  in  comparison 
to  the  much  stronger  no-phonon  paraexciton  line,  comes  from  orthoexcitons 
created  by  Auger  decay  of  paraexcitons  in  the  well,  as  discussed  later. 

The  dependence  of  the  paraexciton  energy  (position  of  the  Xp  line)  on 
applied  stress  is  shown  in  Fig.  4.  The  stress  is  calibrated  using  the 


measured  energy  position  of  the  down-shifted  orthoexciton  line  and  the  data  of 
13 

Waters  et  al.  who  measured  the  orthoexciton  energy  as  a  function  of  applied 
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uniaxial  stress.  The  solid  curve  is  a  fit  of  Eq.  15  to  the  data,  yielding  the 
values  of  A,  B,  and  C  given  in  the  figure  caption.  The  initial  rise  of  the 
paraexciton  energy  is  due  to  the  Bo  term  but  the  Co  dominates  at  high 
stresses  to  give  an  energy  minimum  at  the  SSM. 

In  the  present  experiment,  a  glass  plunger  with  radius  of  curvature 
■  2.5  cm  was  pressed  against  a  (100)  face  of  a  Cu20  crystal  with  dimensions 

3 

1.5  x  1.5  x  1.5  mm  .  Figure  5  is  a  photograph  of  the  stressed  crystal  at  T  - 
1.5  K  with  Ar+  laser  surface  excitation  (X  -  5145A).  The  bright  spot  just 
below  the  top  of  the  crystal  is  luminescence  due  to  excitons  confined  near  the 
SSM.  In  this  case,  the  excitons  must  drift  about  1  mm  from  the  surface  into 
the  potential  well.  The  drift  path  is  barely  visible.  The  luminescence  above 
the  sample  is  due  to  a  reflection  from  the  glass  plunger. 


IV.  Luminescence  Lineshapes  and  Spatial  Profiles 

To  obtain  the  proper  luminescence  lineshape  of  excitons  confined  to 
parabolic  well,  one  must  collect  light  from  all  parts  of  the  well  while 
maintaining  high  spectral  resolution.  In  this  work,  a  method  of  "dynamic 
spatial  integration”  is  used  to  obtain  such  a  luminescence  spectrum.  A  6x 
magnified  image  of  the  sample  is  projected  onto  the  entrance  slit  of  a  1  meter 
spectrometer.  (See  Fig.  1  of  Ref.  15.)  A  1  on  wide  aperture  is  used  to 
select  spatially  the  well  area  in  the  vertical  dimension.  The  size  of  the 
entrance  slit  is  kept  small  in  order  to  have  a  high  spectral  resolution;  but 
this  narrow  slit  could  sample  only  a  fraction  of  the  well.  To  integrate  over 
all  portions  of  the  well,  the  image  of  the  sample  is  scanned  horizontally 
across  the  slit  at  a  5  to  10  Hz  repetition  rate  using  a  triangular  waveform  to 
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drive  Che  galvanometer  so  that  Chere  is  equal  integration  times  for  all 
portions  of  the  well.  The  total  integration  time  is  >  1  second  so  that  many 
passes  of  the  image  are  sampled  at  each  wavelength. 

The  next  data  to  be  described  have  been  obtained  by  pumping  with  a 
tunable  dye  laser  the  orthoexciton  "below”  the  well,  resulting  in  subsequent 
particle  drift  upwards  into  the  well  (see  inset  of  Fig.  6).  The  wavelength  of 
the  dye  laser  is  tuned  to  the  low  energy  edge  of  the  locally  down-shifted 
orthoexciton  phonon-assisted  absorption  band,  resulting  in  an  excited  region 
of  the  crystal  about  0.5  mm  in  length.  Incident  laser  power  is  low,  about  0.3 
mW  CW,  with  about  10Z  absorption.  Each  absorbed  photon  is  converted  into  an 
orthoexciton  and  a  110  cm~*  optical  phonon.  The  orthoexcitons  rapidly  down- 
convert  to  paraexcitons  which  then  drift  upwards  into  the  well. 

In  Fig.  6  no-phonon  paraexciton  lineshapes  taken  at  three  different 
temperatures  are  displayed.  Also  shown  are  fits  of  Eq.  9  to  the  data  assuming 
nondegeneracy  (1(E)  *  E1^2  exp(-E/kgT))  and  using  temperature  as  the  main  fit 
parameter.  (The  two  other  fit  parameters  are  peak  intensity  and  energy 
offset.)  Both  fit  and  bath  temperatures  are  given. 

Phonon-assisted  paraexciton  lineshapes  were  also  taken  under  Identical 
excitation  conditions.  Shown  in  Fig.  7  are  phonon-assisted  lineshapes  taken 
at  the  same  bath  temperatures  as  in  Fig.  6.  These  spectra  are  difficult  to 
obtain  because  of  the  weak  signal  and  the  proximity  of  phonon-assisted 
orthoexciton  lines.  The  applied  stress  for  all  the  data  described  in  this 
section  (o  **  3.6  kbar)  is  chosen  so  that  the  paraexciton  phonon-assisted  line 
is  optimally  separated  from  the  nearby  orthoexciton  phonon-assisted  lines. 

Also  shown  in  Fig.  7  are  fits  of  Eq.  6  to  the  data  assuming  nondegeneracy 
(i.e.,  1(E)  «  E  exp(-E/kgT)).  Again,  temperature  is  the  main  fit  parameter. 
First  consider  the  bottom  spectrum  in  Fig.  7.  The  fit  temperature  of  3.0  K  is 


17 


higher  chan  Che  fic  temperature  of  2.3K  for  Che  no-phonon  line.  There  is  a 
similar  discrepancy  for  Che  middle  and  cop  apecCra  in  Figs.  6  and  7.  Due  Co 
che  ambiguiCy  of  Che  baseline  for  Che  phonon-asaisced  line  ac  higher  tempera- 
Cures,  che  cemperaCure  derived  from  che  Cop  speccrum  in  Fig.  7  has  same 
uncerCainCy  (abouC  ±  1.5  K). 

Clearly,  che  procedure  used  for  fleeing  che  daca  of  Figs.  6  and  7  is 
inadequaCe,  since  che  direcc  and  phonon  assisced  emission  lines  should  yield 
che  same  CemperaCure,  which  is  noc  che  case  here.  What  could  be  causing  a 
difference  between  che  fit  temperatures  of  Che  no-phonon  and  phonon-assisted 
spectra?  One  possibility  is  co  introduce  quantum  degeneracy,  which  will 
affect  Che  no-phonon  lineshape  much  more  than  che  phonon-assisted  line,  as 
explained  in  the  previous  seccion.  Shown  in  Fig.  8  are  the  same  data  as  in 
Fig.  6  but  che  theoretical  fits  to  che  data  (the  solid  dots)  are  obtained  from 
Eq.  9,  assuming  some  degree  of  degeneracy.  The  temperatures  introduced  for 
bath  temperatures  <  3  K  are  fixed  by  Che  phonon-assisted  lineshape  fits  of 
Fig.  7.  The  chemical  potential  is  then  the  main  fitting  parameter  in  Fig. 

8.  The  temperature  and  chemical  potential  for  the  theoretical  fits  are  given 
at  the  left  of  each  spectrum.  The  top  spectrum  is  fit  by  allowing  both  T  and 
p  to  vary.  The  fit  temperature  turns  out  to  be  the  same  as  that  shown  for  the 
top  spectrum  of  Fig.  7.  For  comparison,  nondegenerate  llneshapes  with  the  same 
effective  temperature  as  for  the  degenerate  fits  are  also  shown  in  Fig.  8. 

So  far  the  assumption  of  Bose  degeneracy  would  explain  both  the  no-phonon 
and  phonon-assisted  luminescence  llneshapes.  The  next  question  is  whether  the 
exciton  densities  are  high  enough  to  justify  the  assumption  of  Bose 
degeneracy.  In  order  to  estimate  the  density,  the  relation  N  *  Gt^  is  used, 
where  N  is  the  number  of  excitons,  G  is  the  generation  rate  of  excitons,  and 


Is  Che  Inldal  decay-time  of  the  excltons.  Knowing  the  incident  power  of 
about  0.3  aW  with  10Z  absorption  and  an  initial  decay-time  of  order  0.3  ps, 
then  assuming  100Z  exciton-creation  and  well-collection  efficiency,  one 
obtains  the  estimated  number  of  excltons  in  the  well  to  be  «  2.8  x  10^.  The 

volume  occupied  is  obtained  from  the  measured  spatial  FWHM,  here  about  100 

—7  3 

ua.  This  gives  s  volume  of  -  5  x  10  cm  .  The  estimated  exciton  density  is 
13  —3 

then  •  5.6  x  10  cm  .  This  is  more  than  three  orders  of  magnitude  below  the 
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critical  density  obtained  from  Eq.  4.13,  which  gives  n£  »  10  cm  at  T  -  2 
K.  For  this  reason,  we  are  lead  to  the  conclusion  that  the  existence  of  Bose 
degeneracy  at  this  low  excitation  level  is  highly  unlikely.  Thus,  the  above 
analysis  illustrates  the  danger  in  using  only  spectroscopic  Information  to 
determine  the  degeneracy  of  the  excltonic  system. 

What  could  be  an  alternate  explanation?  A  plausible  explanation  is  as 
follows:  Consider  the  assumptions  made  in  deriving  the  expressions  for  the 
llneshapes  and  spatial  profiles.  The  matrix  element  was  assumed  Independent 
of  spatial  position.  Yet,  it  Is  known  that  the  Intensity  of  the  no-phonon 
paraexciton  line  is  dependent  on  the  applied  stress.  Krelngol'd  and  Makarov^ 

2 

reported  the  Intensity  of  the  no-phonon  paraexciton  line  to  behave  as  o  where 
o  is  the  magnitude  of  the  applied  stress.  The  paraexcltons  in  a  stress 
Induced  potential  well  do  experience  a  range  of  stresses.  The  Intensity  of 
tha  phonon-assisted  line  is  not  greatly  affected  by  stress;  hence,  it  is 
assumed  to  have  a  constant  matrix  element  over  the  range  of  stresses  In  the 
well. 

2 

However,  the  o  intensity  dependence  can  affect  the  no-phonon 
luminescence  llneshape  emitted  from  the  well.  Eq.  7  must  be  modified  to 

Include  the  spatial  dependence  of  the  matrix  element.  Let  the  Intensity  of 

2 

the  paraexciton  be  given  by  Ip  -  So  .  If  the  magnitude  of  the  stress  st  the 
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SSM  Is  oQ,  then  the  stress  distribution  around  the  SSM  can  be  expressed  as 


a  -  oo  - 


-Sr 


(16) 


where  S  is  the  second  derivative  of  stress  with  respect  to  r,  evaluated  at  r  ■ 
0  (the  center  of  the  well).  Expanding  I  about  oQ,  one  obtains 


I 

P 


0o2  +2 $a  (a  -  a  )  +  20(o  -  a  )2 
o  o  o  o 


(17) 


Substituting  Eq.  16  into  Eq.  17,  one  obtains  the  spatial  dependence  of  the 
paraexciton  intensity  (which  is  proportional  to  the  square  of  the  matrix 
element).  Eq.  7  now  becomes 


I(E,r)  « 


I So2  -  20a  Sr2 
o _ o 

«'I(Eloc  - 


*  28S2r4]  SCE 


loc 


-loc’V  -  1 


(18) 


Integrating  over  all  space  now  yields 


1(E) 


E 


1/2 


2nE3/2  +  2n2E5/2 


exp[(E  -  u )/k  T]  -  1 


(19) 


where  n  ■  S/(oQa).  The  second  and  third  terms  in  thj  numerator  represent  the 

correction  to  Eq.  9  due  to  the  stress  dependence  of  the  intensity.  One  can 

estimate  the  parameter  n  knowing  that  S  -  a  /z  2  where  z  is  the  distance  of 

o  o  o 
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the  SSM  from  the  top  surface  (see  Fig.  5).  Using  the  values 

°o  “  3,6  kb*r  *nd  *o  *  0,25  ■“»  009  *et*  S  -  57  kbar/mm2.  With  a  typical 
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2  -1  -1 
a  *  50  meV/nm  ,  this  yields  n  “  0.3  meV  or  300  eV  .  Equation  19  is  shown 

as  the  solid  curve  in  Fig.  9  for  T  -  7.5  K  and  n  "  300  eV~*  and  assuming 

classical  statistics  (u  »  k  T).  The  other  curves  are  the  components  of  this 

B 

lineshape  corresponding  to  the  different  terms  in  Eq.  19  taken  separately  with 

the  occupation  probability  factor,  exp[-E/k  T].  They  show  that  the  stress- 

B 

dependent  Intensity  correction  is  significant.  The  open  circles  are  the  E*^ 

3/2  5/2 

term,  the  open  squares  are  the  E  term,  and  the  crosses  are  the  E  term. 
The  relative  magnitudes  of  the  plots  correspond  to  the  relative  multipliers 
for  each  term.  Comparing  the  solid  curve  with  the  open  circle  curve,  one  sees 
that  Eq.  19  is  narrower  than  Eq.  9  for  a  given  temperature. 

Shown  in  Fig.  10  are  the  no-phonon  lineshapes  with  Eq.  19  fit  to  the 
data.  Both  the  temperature  and  the  parameter  n  are  allowed  to  vary 
lnltally.  It  is  found  that  n  *  300  eV  *  gives  fairly  consistent  results  for 
all  the  bath  temperatures,  with  the  fit  temperature  close  to  the  phonon- 
assisted  lineshapes  at  the  lower  temperatures.  Slight  differences  between  the 
fit  temperature  of  the  no-phonon  and  the  phonon-assisted  lineshapes  still 
exist  for  bath  temperatures  of  2  K  and  below.  Basically,  the  spectral 
inconsistencies  seem  to  be  resolved  by  including  the  stress  dependence  of  the 
no-phonon  luminescence  Intensity.  The  Maxwell-Boltzmann  statistics  used  are 
consistent  with  the  estimated  density. 

Consider,  now  the  spatial  profile  taking  into  account  the  stress- 
dependent  Intensity  of  the  no-phonon  line.  One  must  Integrate  Eq.  18  with 
respect  to  energy.  This  yields 


1  -  2 nor 2  ♦  2n2q2r4 
exp((ar2  -  u )/k^Tj  -  1 


(20) 


Assuming  nondegeneracy,  and  then  integrating  over  the  yz  plane,  one  obtains 


/ 
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k  T  k  T 

I(x)  <*  {1  -  2lna(x2  +  --)]  +  2  [na(x2  +  -^-)]2}  exp(-ax2/k  T)  .  (21) 

a  a  d 

Shown  In  Fig.  11  are  spatial  profiles  of  the  no-phonon  line  taken  at  the  same 

bath  temperatures  as  the  above  described  lineshapes.  The  solid  dots  are  fits 

of  Eq.  21  to  the  data  using  n  “  300  eV-*  and  the  temperatures  from  the  fits  in 

Fig.  10.  The  well  parameter,  a,  was  allowed  to  vary.  It  is  found  that 
2 

a  *  53  meV/mm  gives  the  best  fit  to  spatial  distributions  at  all 
temperatures.  One  sees  that  there  is  excellent  agreement  and  consistency 
among  the  spatial  profiles  and  the  luminescence  lineshapes.  Note  that  for  all 
the  data  described  above,  the  same  values  of  n  and  o  is  used  for  all 
temperatures. 

One  conclusion  that  must  be  drawn  from  this  data  is  that  the  paraexcitons 
in  the  parabolic  well  do  not  seem  to  thermalize  completely  with  respect  to  the 
lattice  within  their  lifetime.  This  is  in  contrast  with  data  taken  on 
crystals  under  zero  applied  stress,  where  thermalized  paraexciton  lineshapes 
are  observed  even  for  much  higher  excitation  conditions.  Notice  that  the  fit 
temperature  in  Fig.  10  are  always  greater  than  the  bath  temperature,  and  the 
disparity  decreases  as  the  bath  temperature  is  lowered.  This  is  opposite  to 
what  one  might  expect  based  on  the  fact  that  the  phonon  scattering  time 
Increases  as  temperature  is  lowered,  which  should  result  in  longer 
thermallzatlon  times  at  lower  temperatures.  There  is,  however,  an  additional 
"thermal izat ion"  that  must  take  place  in  this  case,  i.e.,  a  spatial 
thermallzatlon.  The  paraexcitons  must  fill  the  well  volume,  creating  an 
equilibrium  spatial  distribution  under  the  application  of  this  external 
potential.  At  the  higher  temperatures,  the  equilibrium  size  of  the  spatial 
profile  is  larger  than  at  low  temperatures  (see  Eqs.  2  and  12).  Also,  the 
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diffusivity  of  Che  paraexcitons  at  high  temperatures  Is  much  smaller  than  at 
the  lowe  temperatures*^;  therefore,  it  takes  longer  for  the  excitons  to  fill 
Che  well  at  the  higher  temperatures.  These  factors  would  produce  a  larger 
disparity  between  the  fit  and  bath  temperatures  at  the  higher  bath 
tempe ratures. 


V.  Kinetics  of  Ortho-  and  Paraexcitons 


The  previous  data  was  obtained  at  a  relatively  low  level  of  excitation, 
an  absorbed  power  of  *  0.03  mW.  To  increase  the  density  of  excitons,  one 
naturally  should  increase  the  laser  intensity.  In  so  doing,  we  discovered 
that  the  paraexciton  intensity  (assumed  proportional  to  the  number  of 
particles)  increased  in  a  sublinear  fashion.  Another  puzzling  observation  was 
that  an  orthoexciton  luminescence  signal  emanated  from  the  well  region  even 
when  the  excitation  poi  U  v  as  at  a  considerable  distance  away  from  the  well 
(see  Fig.  3).  The  short  orthoexciton  lifetime  (■  3  ns)  should  preclude 

travel  over  macroscopic  distances.  Also,  thermal  population  of  the 
orthoexciton  level  should  be  negligibly  small:  exp(-AE/kgT)  «  7  x  10-2*  for 
an  ortho-para  splitting  AE  ■>  8  meV  and  temperature  T  «  2  K. 

Shown  in  Fig.  12  is  the  relative  intensity  of  both  ortho-  and  paraexciton 
luminescence  from  the  well  as  a  function  of  absorbed  laser  power,  P.  The 
excitation  configuration  is  the  same  as  that  described  in  the  previous  section 
—  the  laser  absorption  occurs  in  a  region  separated  from  the  well  by  about 
400  pm  and  the  excitons  drift  into  the  well.  The  intensity  of  the 
orthoexciton  was  obtained  from  the  110  cm-*  optical  phonon  replica 
(Xq  -  ri2).  Notice  the  strongly  sublinear  behavior  of  the  paraexciton 
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Intensity  at  high  excitation:  it  is  nearly  a  P  dependence.  The  laser  beam 
was  mechanically  chopped  with  a  4%  duty  factor  at  the  higher  powers  to  avoid 
crystal  heating.  The  power  dependence  of  the  orthoexciton  intensity  is 
significantly  different  from  that  of  the  paraexciton.  In  fact,  the 
orthoexciton  intensity  seems  to  behave  approximately  as  the  square  of  the 
paraexciton  Intensity. 

The  power  dependences  of  the  ortho-  and  paraexciton  intensities  can  be 
modeled  by  rate  equations.  There  are  two  spatial  regions  that  must  be 
considered:  the  region  where  the  laser  light  is  absorbed  and  the  region  of 

potential  well.  First  consider  the  pumped  region  away  from  the  well.  Let  Nq 
and  Np  be  the  number  of  ortho-  and  paraexcitons  in  this  region.  The 
orthoexcitons  are  generated  at  a  rate,  G,  determined  by  the  number  of  photons 
absorbed  per  unit  time,  t.  The  resulting  equations  are 
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where  t  is  the  radiative  lifetime  of  the  paraexciton,  t,  is  the  down- 
P  dc 

conversion  time  for  the  creation  of  paraexcitons  from  orthoexcitons.  The 
aN^  term  is  due  to  an  Auger  decay  process19  in  which  two  paraexcitons  collide 
resulting  in  the  loss  of  one  electron-hole  pair  and  the  ionization  of  one 
paraexciton  into  a  free  electron  and  a  free  hole;  therefore,  the  two 
colliding  paraexcitons  are  lost  in  this  process.  The  free  electron  and  free 
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hole  are  assumed  to  form  an  orthoexclton;  hence,  this  is  a  source  of 
orthoexcitons  (third  term  on  the  right  hand  side  of  Eq.  22b)  but  has  the 
weighting  1/2  in  the  equation  for  the  number  of  orthoexcitons.  It  is  assumed 
that  the  number  of  orthoexcitons  is  small  in  comparison  with  the  number  of 
paraexcitons ,  i.e.,  there  are  no  terms  involving  collisions  between  two 
orthoexcitons  or  between  an  ortho-  and  a  paraexciton.  This  is  a  reasonable 
assumption  in  steady  state  because  the  orthoexclton  down-conversion  rate  is 
much  more  rapid  than  the  paraexciton  decay  rate. 

Now  consider  the  region  of  the  potential  well.  The  equations  for  this 
region  are  virtually  identical  to  those  above  with  the  exception  of  the 
generation  term.  The  source  term  of  paraexcitons  is  assumed  to  be  from 
drifting  into  the  well,  at  a  rate  Gp,  which  is  assumed  proportional  to  Np. 

Let  nQ  and  np  be  the  number  of  ortho-  and  paraexcitons  in  the  well.  The  rate 
equations  for  the  well  region  are  then 
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Again,  the  number  of  orthoexcitons  is  assumed  small  compared  to  the  number  of 
paraexcitons  in  steady  state.  Solving  Eqs.  22  and  23  for  the  number  of 
paraexcitons  in  steady  state  yields 


N  =  — - —  [-1  +  A  +  2aGx2'  ] 
p  axp  p 
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In  steady  state, 
cNp,  where  c  can 


2 

Eq.  23b  shows  that  nQ  is  proportional  to  np  .  Assuming  Gp  = 
be  thought  of  as  the  well  capture  rate,  one  obtains 
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The  solid  curve  marked  1  in  Fig.  12  is  Eq.  25a  with  =  0.5  ps  ,  a  = 

0.05  sec-*,  and  c  =  1  x  10^  sec  _*.  The  value  of  t  is  determined  by  a  low 

P 

power  decay-time  measurement  to  be  described  below.  The  parameters  a  and  c 
are  both  allowed  to  vary  initially.  It  is  found  that  the  shape  of  the  curve 
is  more  sensitive  to  the  value  of  a  than  that  of  c. 

One  can  independently  estimate  the  value  of  c  in  the  following  manner. 

The  product  cNp  is  the  number  of  paraexcitons  per  unit  time  flowing  into  the 
well.  The  number  which  reach  the  well  is  of  order 

Npexp(-ttransit./ip),  with  ttranslc  equal  to  the  amount  of  time  a  particle 
takes  to  reach  the  well  from  the  pumped  region  of  the  crystal.  This  gives  c  = 
(exp  (-ttranslt/Tp  ))/ttransit.  Now  ttranslt  =  Ax/vd,  where  Ax, 
the  distance  traveled,  is  approximately  0.4  mm,  and  v^,  the  drift  velocity,  is 
approximately  1  x  10"*  cm/s.  This  yields  c  “  1  x  10^  sec  -*.  The  parameter  c 
is  fixed  at  this  value  and  the  parameter  a  is  varied  to  give  a  reasonable  fit 
to  the  data.  The  line  marked  2  is  Eq.  25b  with  the  adjustment  of  an  overall 
multiplier  to  give  the  appropriate  vertical  scale. 

The  nearly  behavior  of  the  paraexciton  intensity  in  the  data 

described  above  is  due  to  an  Auger  decay  process  at  two  different  regions  of 


the  crystal.  If  one  were  to  create  the  excitons  directly  in  the  well  area, 
one  region  is  eliminated  and  the  power  dependence  of  the  paraexciton  intensity 
should  be  P*^.  This  is  seen  by  assuming  Gp  «  P  in  Eq.  24b.  Shown  in  Fig.  13 
is  the  extracted  density  of  paraexcitons  in  the  well  when  the  laser  is  pumping 
the  well  area  directly.  This  data  shows  the  predicted  P  '  dependence.  (The 
corresponding  orthoexciton  intensity  data  was  not  taken  at  that  time.)  Thus, 
the  data  from  this  simpler  pumping  configuration  support  the  basic  rate 
equations  we  have  used. 

Another  test  of  this  model  is  provided  by  the  transient  decay  of  the  two 
types  of  excitons.  Shown  in  Fig.  14  are  time  decays  of  the  ortho-  and 
paraexciton  signals  after  resonant  creation  of  orthoexcitons  at  the  "bottom" 
of  the  well  (i.e.,  the  laser  was  tuned  to  the  orthoexciton  quadrupole 
resonance  at  the  SSM).  The  dye  laser  was  cavity-dumped.  After  an  initial 
fast  decay,  the  orthoexciton  signal  decays  at  twice  the  rate  of  the 
paraexciton  signal.  This  is  consistent  with  the  Auger  decay  mechanism  as  can 
be  seen  by  substituting  into  Eqs.  22  a  solution  of  the  form  N  =  N°exp 
(-t/Te^)  for  both  the  ortho-  and  paraexciton  number  (a  different  for 

each).  The  generation  rate  G  is  zero  at  times  after  t  =  0.  With  the 
additional  information  that  x^c  is  of  order  a  few  nanoseconds,^  and  that  the 
measured  decay-times  are  a  few  tens  to  hundreds  of  nanoseconds,  t(jc  can  then 
be  ignored  in  comparison  to  *eff  for  the  orthoexcitons.  This  easily  yields 
the  results  that  the  decay-time  for  the  paraexcitons  should  be  twice  that  of 
the  orthoexcitons. 

The  inset  to  Fig  14  shows  the  measured  paraexciton  decay-times  for  three 
incident  powers.  Estimation  of  absorbed  power  is  impossible  because  of  the 
pumping  configuration.  Both  the  initial  and  the  later  decay-times  are 
given.  The  decay-times  decrease  as  the  power  is  lowered,  as  expected  for  an 
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Auger  loss  mechanism.  At  the  lowest  power,  there  is  little  difference  between 

the  fast  and  the  slow  decay-times.  It  is  this  low  power  decay  curve  which  is 

the  source  of  the  value  of  t  used  above. 

P 


VI.  Thermalization  of  Excitons  in  a  Parabolic  Well 


A  basic  question  that  is  raised  by  all  the  above  data  is  whether 
paraexcitons  energetically  and  spatially  thermalize  in  the  parabolic  well.  To 
gain  more  insight  into  what  is  happening,  time-resolved  spectra  and  spatial 
profiles  were  taken  under  various  excitation  conditions.  Before  discussing 
data  taken  from  the  well  area,  contrasting  data  taken  under  high-powered 
pulsed  excitation  with  no  stress  applied  to  the  crystal  will  be  discussed 
briefly. 

Shown  in  Fig.  15  are  some  time-resolved  spectra  of  the  phonon-assisted 
paraexclton  line  following  pulsed  excitation  of  the  crystal  at  a  bath 
temperature  of  2  K.  The  peak  incident  power  was  about  3  watts  and  the  pulse 
length  was  about  15  ns  (cavity-dumped  dye  laser).  The  laser  was  tuned  into 
the  phonon-assisted  orthoexciton  absorption  band  with  an  absorption  length  of 

about  0.3  mm.  The  time  t=0  corresponds  to  the  peak  of  the  laser  pulse. 

1/2 

Theoretical  fits  to  the  spectra  are  an  E  '  exp(-E/kgT)  Maxwell-Boltzmann 
distribution  with  temperature  as  a  variable  parameter.  The  fit  temperatures 
are  shown  at  the  left.  Note  that  the  paraexcitons  display  a  nearly 
thermalized  distribution  of  particle  energies  at  short  delay  times  after  the 
intense  laser  pulse. 

Now  consider  data  taken  from  the  well  under  rather  mild  excitation 
conditions.  The  laser  was  tuned  to  the  orthoexciton  quadrupole  resonance  at 
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Che  3SM  ("bottom"  of  the  well).  The  incident  peak  power  was  about 
0.5  W.  Estimation  of  the  absorbed  power  was  impossible  because  of  the  pumping 
configuration,  but  it  should  be  a  very  small  fraction  of  the  incident  power 
due  to  the  weakness  of  the  quadrupole  transition  and  the  short  distance  over 
which  the  absorption  takes  place.  Shown  in  Fig.  16a  are  the  resulting  time- 
resolved  spectra.  The  time  delay  between  displayed  spectra  is  80  ns.  At 
early  times  the  spectra  are  narrow  due  to  the  narrowness  of  the  spatial 
distribution.  As  the  excitons  diffuse  outward  to  fill  the  well,  the  spectra 
broaden  correspondingly.  Shown  in  Fig.  16b  is  an  effective  temperature,  which 
is  a  measure  of  the  width  of  the  spectrum,  as  a  function  of  time.  This 
effective  temperature  was  obtained  by  assuming  a  simple  E  exp(-E/kgTe^j ) 
lineshape  which  yields  l.8kgTg££  for  the  value  of  the  FWHM.  No  attempts  were 
made  to  correct  for  the  stress  dependent  intensity  of  the  paraexciton  no¬ 
phonon  transition.  The  feature  of  note  is  that  at  later  times  the  excitons 
seem  to  exhibit  an  effective  temperature  higher  than  that  of  the  bath,  and  the 
thermalization  seems  to  be  very  slow.  This  is  in  direct  contrast  to  the  more 
rapid  thermalization  observed  at  zero  stress  (Fig.  15). 

Shown  in  Fig.  17a  are  time-resolved  spatial  profiles  taken  under 
virtually  identical  experimental  conditions  as  above.  Again,  the  time  delay 
between  displayed  profiles  is  80  ns.  Initially  the  spatial  distribution  has  a 
narrow  width  as  determined  by  the  size  of  the  laser  beam.  As  the  excitons 
diffuse  outward  to  fill  the  well,  the  width  grows.  Shown  in  Fig.  17b  is  the 
square  of  the  FWHM,  A,  as  function  of  time.  The  width  of  the  spatial  profile 
does  decrease  slightly  at  longer  times,  possibly  indicating  a  cooling  process. 

One  can  model  this  diffusive-like  expansion  process  by  assuming  that  the 
mean  free  path  of  excitons  between  collisions  with  phonons  is  small  compared 
to  the  macroscopic  dimension  of  the  potential  well.  This  allows  one  to  model 


the  external  potential  as  providing  a  spatially  varying  motive  force  on  the 
excitons  which  determines  a  "local''  drift  velocity.  Along  one  dimension,  an 
external  potential  of  ax  gives  a  force,  F(x)  ■  -2ax.  This  results  in  a 
"local"  drift  velocity,  v(x)  -  F(x)  x/m*,  where  T  is  the  phonon  scattering 
time.  Particle  conservation,  with  a  sink  due  to  the  finite  lifetime  of  the 
excitons,  then  yields  the  following  differential  equation 
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where  y  *  2at/m*,  n  ■  n(x,t)  is  the  local  "  ensity”  of  particles,  D  is  the 
diffusion  constant,  and  is  the  particle  lifetime.  The  solution  to  this 
equation  which  is  a  6-function  at  x  *  0,  t  =  0  is 
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The  FWHM  of  this  distribution  is  given  by 

A2(t)  -  2.77(4)  -2_  (i  -  e~2Yt).  (28) 

The  solid  line  in  Fig.  17b  is  a  fit  of  Eq.  28  to  the  data  allowing  D, 

y,  and  the  "zero”  of  time  to  vary.  (This  is  required  because  the  spatial 

distribution  has  a  finite  size  when  created  by  the  laser.)  The  resultant 

2 

parameters  are  given  in  the  figure  caption.  A  value  of  D  -  190  cm  / s  is 
obtained  from  the  fit.  This  reasonably  consistent  with  the  measured^  drift 


mobility  of  paraexcltons  at  a  lattice  temperature  of  2  K  and  at  a  high  applied 
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stress.  Specifically,  u  -  7.5  x  10^  cm^/eV-s  Implies  D  ■  130  cm^/s  from  D 

»  ukgT. 

Finally,  experiments  have  been  performed  at  the  highest  paraexciton 
densities  possible  with  the  apparatus.  The  laser  wavelength  and  spatial 
position  are  adjusted  to  give  the  maximum  no-phonon  paraexciton  signal  from 
the  well.  The  optimal  excitation  configuration  is  very  similar  to  that  used 
to  obtain  the  potential-well  luminescence  lineshapes  described  in  section 
IV.  The  differences  are  that  the  laser  wavelength  is  shorter,  resulting  in  an 
absorption  region  nearly  the  width  of  the  sample,  and  the  vertical  position  of 
the  pump  beam  is  closer  to  the  well  area.  Maximum  available  cavity-dumped  dye 
laser  power  is  used,  about  5  W  peak  incident  on  the  sample.  The  fraction  of 
absorbed  power  is  in  the  range  10-20%,  as  determined  from  transmission 
measurements  in  similar  configurations. 

Shown  in  Fig.  18  are  time-resolved  spectra  taken  from  the  well  area  under 
these  extreme  excitation  conditions.  The  time  delay  between  displayed  spectra 
is  80  ns.  The  very  strange  shape  of  the  spectra  at  early  times  (bottom 
spectrum)  is  due  to  paraexcitons  decaying  on  their  way  into  the  well.  Since 
they  decay  in  a  region  of  the  crystal  at  a  lower  stress,  the  emitted  light 
will  be  of  higher  energy.  Note  the  long  time  (■  200  ns)  it  takes  the  excitons 
to  finally  arrive  in  the  well  and  display  a  Maxwellian  lineshape.  This  is 
Interpreted  as  due  to  a  true  lattice  heating  at  these  high  powers.  Similar 
data  at  a  factor  of  30  lower  power  shows  a  much  more  rapid  arrival  in  the 
well,  thus,  displaying  a  Maxwellian  lineshape  at  earlier  times.  There  still 
seems  to  be  a  slow  therraalization  taking  place.  The  effective  temperature  of 
the  top  spectrum  in  Fig.  18  as  determined  from  its  FWHM  is  about  3.5  K  at  a 
bath  temperature  of  1.3  K. 

The  data  described  in  this  paper  show  a  rather  slow  thermalizat ion  of 


paraexcltons  in  a  stress-induced  parabolic  potential  well.  The  time-resolved 
spatial  profiles  indicate  that  the  excltons  are  able  to  "fill"  the  well  in 
times  <  200  ns,  due  to  their  rather  large  diffusion  constant.  But  the 
effective  temperature  derived  from  the  spectral  distribution  remains  1  to  3  K 
above  the  lattice  temperature  over  the  lifetime  of  the  particle,  with  the 
amount  of  heating  dependent  on  excitation  power.  In  contrast,  when  no  stress 
is  applied  to  the  crystal,  the  paraexcltons  thermalize  much  more  quickly,  even 
under  fairly  higher  laser  excitation.  What  is  the  major  difference  between 
these  two  situations? 

It  Is  known  that  external  uniaxial  stress  causes  the  energy  separation 
between  the  paraexciton  and  the  down-shifted  component  of  the  orthoexciton  to 
decrease,  and  paraexcltons  are  produced  by  the  down  conversion  of 

orthoexcitons.  Ortho-para  down  conversion  has  been  studied  to  some 

7  1 6  2C 

extent  ’  ’  ,  but  there  is  still  no  consensus  as  to  the  exact  processes 

involved.  With  no  stress  applied,  96  cm  *  (12  tneV)  of  energy  must  be  given 

off  during  the  down-conversion.  This  may  occur  via  the  emission  of  an  87  cm  ^ 

-  20 

optical  phonon  of  symmetry  and  an  acoustic  phonon.  However,  the 
predicted  temperature  dependence  for  the  rate  of  this  process  was  not  observed 
by  Weiner  et  al. ^ 

The  participation  of  this  87  cm  ^  optical  phonon  may  still  provide  a  key 
to  explaining  the  difference  in  stress  and  unstressed  thermalization  times. 
When  a  stress  of  only  about  1  kbar  is  applied  to  the  crystal,  the  energy 
separation  between  the  paraexciton  and  the  down-shifted  orthoexciton  levels 
becomes  less  than  87  cm  \  which  could  alter  the  down-conversion  process, 
possibly  resulting  in  a  hotter  initial  distribution  of  paraexcltons  which  must 
then  thermalize  down  to  the  lattice  temperature. 
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VII.  Summary  and  Conclusions 

Bxcitons  In  a  direct-gap  semiconductor  have  been  successfully  strain- 
confined  for  the  first  time.  In  contrast  to  the  cases  of  the  indirect  gap 
materials  Si  and  Ge,  where  the  potential  well  arises  from  a  splitting  of  a 
degenerate  conduction  band,  here  the  stress-induced  potential  well  at  the  SSM 
is  due  to  an  interaction  between  the  valence  bands  causing  the  symmetry¬ 
lowering  component  of  the  applied  stress  to  be  the  dominating  factor. 

Specific  predictions  are  made  for  the  luminescence  lineshapes  and  spatial 
profiles  of  luminescence  emitted  by  strain-confined  excltons.  It  is  found 
that  as  the  system  of  excltons  is  allowed  to  become  quantum-degenerate,  the 
phonon-assisted  line  gives  little  indication  of  the  occurence  of  the  Bose 
degeneracy;  however,  the  no-phonon  line  should  show  significant  changes  as 
degenerate  Bose  statistics  sets  in.  For  this  reason,  the  existence  of  a  no¬ 
phonon  transition,  such  as  is  the  case  in  a  direct-gap  semiconductor,  may 
prove  quite  beneficial  when  trying  to  observe  quantum  statistical  behavior  for 
excltons  in  a  parabolic  well. 

Spectral  lineshapes  of  paraexclton  luminescence  emitted  from  the  well 
area  for  both  the  no-phonon  and  the  phonon-assisted  lines  of  C^O  have  been 
studied.  Discrepancies  between  fit  temperatures  of  the  no-phonon  and  the 
phonon-assisted  lines  are  found  under  the  assumption  of  a  simple  lineshape 
theory  and  Maxwell-Boltzmann  statistics  for  the  decaying  excitons.  On  the 
other  hand,  both  the  direct  and  phonon-assisted  paraexclton  luminescence 
lineshapes  could  be  accounted  for  with  the  same  effective  temperature  in  the 
simple  llne8hape  theory  by  Invoking  Bose  degeneracy  of  the  decaying 
paraexcltons.  However,  the  estimated  densities  were  much  too  low. 

Invalidating  this  assumption.  A  plausible  alternate  explanation  incorporates 
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the  stress  dependence  of  the  paraexclton  no-phonon  transition  and  yields 
consistent  results  in  the  classical  statistical  limit.  There  was  a  good 
correspondence  between  the  measured  no-phonon  lineshapes  and  spatial  profiles. 

The  intensity  of  ortho-  and  paraexclton  luminescence  emitted  from  the 
well  region  was  also  studied  as  a  function  of  laser  power.  The  paraexclton 
intensity  displayed  a  very  subllnear  dependence  at  high  laser  powers.  This 
sublinear  behavior  was  attributed  to  an  Auger  decay  mechanism  resulting  from 
collisions  among  the  paraexcltons.  Due  to  the  thermalization  of  the  free 
electron  and  hole  created  in  the  Auger  decay  process,  there  was  an  additional 
source  of  orthoexcltons ;  hence,  the  appearance  of  orthoexciton  signal  in  the 
well  even  if  the  excitation  source  was  Located  at  a  large  distance  from  the 
we  1 1 . 

Lastly,  time-resolved  data  were  presented  which  explored  the 
thermalization  of  paraexcltons  in  the  parabolic  well  under  various  excitation 
conditions.  Even  under  fairly  mild  excitation  conditions  there  seemed  to  be  a 
very  slow  thermalization  taking  place.  This  was  contrary  to  the  more  rapid 
thermalization  observed  for  the  paraexcltons  in  crystals  with  no  stress 
applied.  Also,  under  Intense  laser  excitation,  significant  lattice  heating 
occurred . 

In  attempting  to  achieve  Bose-Einstein  condensation  of  long-lived 
paraexcltons  in  (^O,  we  have  thus  encountered  two  barriers.  The  principal 
one  is  the  limitation  of  exciton  density  due  to  Auger  collisions  even  at  quite 

*  i  n 

small  particle  densities  (>  101  cm"J).  A  secondary  barrier  is  the  rather 
slow  thermalization  of  excitons  in  the  well.  These  formidable  difficulties 
must  be  surmounted  before  Bose-Einstein  condensation  of  long-lived 
paraexcltons  confined  to  a  strain  well  can  be  realized  in  C^O. 
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Fig.  1 


Fig.  2 


Fig.  3 


Figure  Captions 

Theoretical  predictions  for  the  luminescence  lineshape  (part  a)  and 
the  spatial  profile  (part  b)  of  a  phonon-assisted  emission  line 
from  excitons  confined  to  a  parabolic  potential  well.  The 
different  curves  are  for  varying  degrees  of  quantum  degeneracy. 

Note  the  little  change  that  occurs  in  both  the  luminescence  line- 
shape  and  the  spatial  profile  when  some  quantum  degeneracy  is 
assumed.  The  bottom  trace  in  each  part  schematically  includes 
macroscopic  occupancy  of  the  ground  state  when  p  =  0. 

Theoretical  predictions  for  the  luminescence  lineshape  (part  a)  and 

spatial  profiles  (part  b)  of  a  no-phonon  emission  line  from 
excitons  confined  to  a  parabolic  potential  well.  The  different 

curves  are  for  varying  degrees  of  quantum  degeneracy.  There  are 
significant  changes  in  both  the  luminescence  lineshape  and  the 
spatial  profile  when  some  degeneracy  is  assumed.  The  bottom  trace 
in  each  part  schematically  includes  macroscopic  occupancy  of  the 
ground  state  when  p  =  0.  Both  these  degenerate  traces  have 
singularities  which  are  shown  as  the  arrows  pointing  upwards. 
Luminescence  spectra  taken  with  no  stress  applied  to  the  crystal 
(lower  trace),  and  with  a  stress  of  about  3.6  kbar  at  the  shear- 
stress  maximum  (upper  trace).  Note  the  different  wavelength  scales 
for  the  two  spectra.  XQ  is  the  no-phonon  orthoexciton  emission  and 
the  other  lines  are  phonon  replicas  of  the  orthoexciton  labeled  by 
the  symmetry  of  the  optical  phonon  emitted,  with  the  exception  of 

the  Xp  -  line  which  is  the  phonon-assisted  paraexciton 
emission.  A  new  line  appears  when  uniaxial  stress  is  applied  to 
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Fig.  4 


Fig.  5 


Fig.  6 


the  crystal  —  the  no-phonon  paraexciton  line  labeled  X^.  This 
line  is  cut  off  at  about  1/8  its  maximum  intensity  and  the  XQ 
-  line  is  cut  off  at  about  1/3  its  maximum  intensity. 

Energy  position  of  the  paraexciton  as  a  function  of  applied  stress, 
o.  The  data  were  obtained  by  using  the  measured  energy  shift  of 
the  orthoexciton  line  as  the  stress  calibrator.  The  fit  to  the 
orthoexciton  versus  stress  data  given  in  Ref.  13  was  used  as  the 
calibrating  curve.  The  solid  curve  is  a  fit  of  Eq.  15  to  the  data 

with  the  values  A  =  16298  cm  * ,  B  =  4.96  cm  ^kbar  *  and  C  =  7.36 
-1  -2 

cm  kbar  .  The  contribution  of  the  separate  terms  of  Eq .  15  are 

2 

shown  as  the  dashed  curves.  The  dominance  of  the  Co  term  at  high 
stresses  makes  the  shear-stress  maximum  an  energy  minimum  for  the 
paraexcitons. 

Photograph  of  a  C^O  crystal  with  stress  applied  by  a  glass 
plunger.  An  argon-ion  laser  excites  the  surface  of  the  sample. 

The  created  excitons  drift  in  the  stress  gradient  to  the  shear- 
stress  maximum  where  they  are  confined.  This  region  is  the  bright 
spot  near  the  top  of  the  sample.  The  luminescence  seen  above  the 
sample  is  due  to  a  reflection  from  the  glass  plunger.  The  sample 

3 

dimensions  are  1.5  x  1.5  x  1.5  mm  . 

No-phonon  paraexciton  luminescence  lineshapes  for  three  bath 

temperatures  taken  using  "dynamic  spatial  integration"  of 
luminescence  emitted  from  the  potential  well  as  explained  in  the 

text.  The  fits  to  the  data  (solid  dots)  are  the  nondegenerate 
limit  of  Eq.  9  (i.e.,  exp(-E/kgT))  convolved  with  the 

resolution  function  shown  at  the  left  of  each  trace.  The  resultant 


fit  temperatures  are  given  at  the  left  and  the  bath  temperatures 
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Fig.  7 


Fig.  8 


Fig.  9 


are  given  at  the  right.  The  magnitude  of  the  applied  stress  was 
about  3.6  kbar. 

Phonon-assisted  paraexciton  luminescence  lineshapes  for  the  same 
bath  temperatures  as  Fig.  6  taken  using  "dynamic  spatial 
integration"  of  luminescence  emitted  from  the  potential  well.  The 
fits  to  the  data  (solid  dots)  are  the  nondegenerate  limit  of  Eq.  6 
(i.e.,  E  exp(-E/kgT))  convolved  with  the  resolution  function  shown 
in  the  upper  right-hand  corner.  Fit  temperatures  are  shown  at  the 
left.  The  "fit"  temperature  of  7.5  K  for  the  top  trace  is  a  bit 
uncertain  due  to  ambiguous  baseline  determination  which  is  a  result 
of  nearby  orthoexcitation  phonon-assisted  lines.  The  determination 
of  this  temperature  is  explained  in  the  text. 

No-phonon  paraexciton  luminescence  lineshapes  again,  the  same  data 
as  in  Fig.  6.  These  fits  to  the  data  (solid  dots)  are  Eq.  9 
assuming  some  quantum  degeneracy.  The  fit  temperatures  for  the 
lower  two  spectra  were  fixed  at  the  values  obtained  for  the  phonon- 
assisted  lineshapes  (see  Fig.  7)  and  the  chemical  potential,  q,  was 
allowed  to  vary.  The  resulting  values  of  p  are  given  at  the  left 
as  a  fraction  of  k^T  where  T  is  the  fit  temperature.  For  the  upper 
spectrum,  both  T  and  p  were  varied  for  the  best  fit.  Shown  as  the 
x's  is  the  corresponding  nondegenerate  lineshape  with  the  same 
temperature  as  the  fit  temperature  given.  The  assumption  here  of 
quantum  degeneracy  is  actually  invalid  as  explained  in  the  text. 
Luminescence  lineshape  for  the  no-phonon  paraexciton  line  taking 
into  account  the  variation  in  emission  intensity  with  stress. 
Maxwell-Boltzmann  statistics  (p  »  k_T)  is  assumed.  The  broken 
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curves  (circles,  boxes  and  x's)  are  the  components  of  the  lineshape 
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Fig. 


Fig. 


Fig. 


Fig. 


as  given  by  the  different  terms  in  Eq.  19,  with  an  inverted  curve 
indicating  a  negative  contribution.  The  resultant  lineshape  is 
shown  as  the  solid  curve.  A  temperature  of  7.5  K  was  used  with  the 
parameter  ri  =  300  eV-^.  Note  that  this  lineshape  is  narrower  at  a 
given  temperature  than  the  lineshape  given  in  Eq.  9. 

10  No-phonon  paraexclton  luminescence  lineshapes  again,  the  same  data 
as  in  Fig.  6.  The  fits  to  the  data  (solid  dots)  are  Eq.  19 
convolved  with  the  resolution  function  shown  at  the  left  of  each 
curve  and  p  »  kgT.  A  single  value  of  n  =  300  eV-*  was  used  for 
all  fits.  Temperature  was  allowed  to  vary  somewhat  but  good  fits 
were  obtained  with  temperatures  the  same  as  or  very  close  to  the 
fit  temperatures  obtained  from  the  phonon-assisted  lineshapes  (see 
Fig.  7). 

11  Spatial  profiles  of  the  no-phonon  paraexciton  line  at  the  same  bath 
temperatures  as  Fig.  10.  The  fits  to  the  data  (solid  dots)  are  Eq. 
21  with  the  same  value  n  =  300  eV~^  and  the  same  fit  temperatures 
as  those  in  Fig.  10.  Shown  in  the  upper  right-hand  comer  is  the 
estimated  spatial  resolution. 

12  Intensity  of  the  no-phonon  paraexciton  line  and  the  phonon 
replica  of  the  orthoexciton  as  a  function  of  absorbed  laser  power. 
The  excitation  configuration  is  shown  in  the  inset.  The  laser  is 
absorbed  in  a  region  below  the  potential  well  and  the  excitons 
drift  upwards  into  the  well  (shown  as  the  dot  near  the  top  of  the 
sample).  The  solid  curves  are  fits  of  Eq.  25  to  the  data. 

13  Estimated  density  of  paraexcitons  as  a  function  of  absorbed  laser 

power  when  exciting  the  well  area  directly.  The  dependence  is 

consistent  with  the  Auger  decay  mechanism  proposed  in  the  text. 
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Fig. 


Fig. 


Fig. 


14  Intensity  of  ortho-  and  paraexciton  luminescence  as  a  function  of 
time  after  resonant  creation  of  orthoexcitons  at  the  shear-stress 
maximum  with  a  cavity-dumped  dye  laser  pulse.  The  decays  are 
nonexponential.  After  an  initial  fast  decay  the  orthoexciton 
signal  decays  at  twice  the  rate  of  the  paraexciton  signal.  The 
lines  through  the  data  are  labeled  with  the  corresponding  1/e 
times.  Shown  in  the  inset  are  the  fast  (solid  squares)  and  the 
slow  (dots)  decay  times  of  the  paraexciton  signal  as  a  function  of 
incident  peak  laser  power.  Note  that  as  the  power  is  lowered,  the 
decay  seems  to  become  much  more  exponential,  i.e.,  one  decay  time. 

15  Time-resolved  spectra  of  the  phonon-assisted  paraexciton  line  with 
no  stress  applied  to  the  crystal  and  intense  cavity-dumped  dye 
laser,  near-surface  excitation.  The  fits  to  the  data  (solid  dots) 
are  Maxwell-Boltzmann  distribution  E1/,2exp(-E/kgT)  convolved  with 
the  resolution  function  shown  in  the  upper  right-hand  corner.  Fit 
temperatures  are  shown  at  the  left  with  the  delay  time,  t,  given  at 
the  right  (t  =  0  defined  as  near  the  center  of  the  laser  pulse). 
Note  that  the  excitons  display  a  near  thermal  distribution  of 
particle  energies  at  short  delay  times  after  the  laser  pulse. 

16  a)  Time-resolved  spectra  of  the  no-phonon  paraexciton  line  after 
resonant  creation  of  orthoexcitons  at  the  shear-stress  maximum. 
There  is  80  ns  delay  between  displayed  spectra,  with  the  bottom 
spectrum  at  the  earliest  time. 

b)  An  effective  temperature  extracted  from  the  full-width  at  half¬ 
maximum  of  the  spectra  in  (a)  assuming  a  standard  E^2exp(- 


E/kgTgf f )  form  for  the  lineshape.  As  the  paraexcitons  expand  to 
fill  the  well,  the  luminescence  lineshape  broadens  correspondingly. 


42 


t  »  0  Is  near  the  peak  of  the  laser  pulse. 

Fig.  17  a)  Time-resolved  spatial  profiles  of  no-phonon  paraexciton 

luminescence  intensity,  I(x),  after  resonant  creation  of  ortho- 
excitons  at  the  shear-stress  maximum.  There  is  80  ns  delay  between 
displayed  profiles. 

2 

b)  Square  of  the  full-width  at  half-maximum,  A  ,  as  a  function  of 

time.  The  solid  line  is  a  fit  of  Eq.  28  to  the  data  with  the  para- 
2  6  “1 

meters  D  =  190  cm  /s,  y  =  5.6  x  10  sec  ,  and  a  "zero"  of  time 
about  7  ns  before  the  peak  of  the  laser  pulse  (see  text). 

Fig.  18  Time-resolved  spectra  of  the  no-phonon  paraexciton  line  emitted 

from  the  well  area  while  pumping  horizontally  below  the  potential 
well  with  fairly  intense  laser  excitation.  The  excitons  drift  into 
the  well.  There  is  80  ns  time  delay  between  displayed  spectra  with 
the  earliest  time  at  the  bottom.  The  strange  shape  of  the  spectra 
at  early  times  is  due  to  the  decay  of  particles  on  their  way  into 
the  well.  Even  for  the  latest  time  displayed,  the  spectrum 
displays  a  distribution  much  hotter  than  the  bath  temperature. 
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